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^ , Abstract. This is the fourth in a series of papers that aim both to provide reasonable orbits for a number of 

' eclipsing binaries and to evaluate the expected performance of Gaia of these objects and the accuracy that is 

I achievable in the determination of such fundamental stellar parameters as mass and radius. In this paper, we 

attempt to derive the orbits and physical parameters for three eclipsing binaries in the mid-F to mid-G spectral 
' range. As for previous papers, only the Hp, Vt, Bt photometry from the Hipparcos/Tycho mission and ground- 

tr^ ^ based radial velocities from spectroscopy in the region 8480-8740 A are used in the analyses. These data sets 

lO ■ simulate the photometric and spectroscopic data that are expected to be obtained by Gaia, the approved ESA 

' Cornerstone mission to be launched in 2011. The systems targeted in this paper are SV Cam, BS Dra and HP Dra. 

SV Cam and BS Dra have been studied previously, allowing comparisons of the derived parameters with those 
^ from full scale and devoted ground-based investigations. HP Dra has no published orbital solution. SV Cam has a 

(3 Lyrae type light curve and the others have Algol-like light curves. SV Cam has the complication of light curve 
anomalies, usually attributed to spots; BS Dra has non-solar metallicity, and HP Dra appears to have a small 
eccentricity and a sizeable time derivative in the argument of the periastron. Thus all three provide interesting 
and different test cases. 
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>y ,. 1. Introduction Straizys (1999), Bienayme & Turon (2002), Vansevicius et 

al. (2002), Munari (2003) and Turon & Ferryman (2004). 

Gaia is a Cornerstone mission re-approved by ESA in May 

2002. It is intended to perform three important tasks: (a) In Paper I of this series, Munari et al. (2001), we be- 
micro-arcsec astrometry, (b) photometry in 4 broad and gan the process of evaluating the expected performance 
11 intermediate passbands and (c) R 11500 resolution by Gaia in dealing with eclipsing binaries and the accu- 
spectroscopy in the AA8480— 8740 A wavelength region, racy with which fundamental stellar parameters such as 
The completeness limit for both astrometry and photom- masses, radii and temperatures can be determined. The 
etry is anticipated to be ~ 20 mag. Each target star number of eclipsing binaries to be discovered by Gaia is ex- 
is expected to be measured around a hundred times dur- pected to be ^^10^, of which 10^ are likely to be double- 
ing the five year mission life-time, in an operational mode lined. This is orders of magnitude larger than all the SB2 
similar to that of Hipparcos. The astrophysical goals and eclipsing binaries so far investigated from ground-based 
technical specifications of the mission can be found, among observations (e.g. Andersen 1991, 2002), with the addition 
others, in ESA's Concept and Technology Study (ESA of stunningly precise astrometric parallaxes to be used in 
SP-2000-4), Gilmore et al. (1998), Ferryman et al. (2001) close-loop iterative refinement of parameters involved in 
and in the proceedings of Gaia conferences edited by the orbital modeling. Furthermore, Gala's discoveries will 

be a boon to ground-based observatories by providing op- 
portunities to follow up discoveries with synoptic, multi- 
Send offprint requests to: U. Munari wavelength campaigns for many decades thereafter. It is 
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Table 1. Eclipsing binary targets. Data from the Hipparcos Catalogue. Bt and Vt are Tycho mean and Hp are 
median magnitude values. 



Names 



Sp. Hp Bt 



Vt 



ajl991.25 

(h m s) 



5,71991.25 



parallax /i^ 
(mas) (mas yr~^) 



Us 

(mas yr~^) 



SV Cam HIP 32015 G5 9.4377 10.204 9.411 06 41 18.89 +82 16 03.8 11.77±1.07 +41.58±0.95 -152.91±1.17 
BS Dra HIP 98118 F5 9.2320 9.638 9.183 19 56 28.79 +73 36 57.6 4.80+0.74 -7.95+0.68 -5.06+0.72 
HP Dra HIP 92835 G5 8.0888 8.679 8.032 18 54 53.46 +51 18 29.1 12.45±0.72 +23.35+0.77 +83.40+0.90 




0500 



8600 



8700 



Fig. 1. Sample of the spectra of each of the variables discussed here. Note the lack of visible features of a secondary 
component in the SV Cam spectrum contrary to some controversial claims in literature. The spectra are shifted to 
heliocentric corrected wavelengths. 



therefore of great interest to test Gaia's capabilities for 
eclipsing binary investigations, with the aim of contribut- 
ing to the fine tuning of mission planning and preparing 
analysis strategies for the massive data flow that will fol- 
low. 

Paper I outlined the framework of the project and 
adopted methodologies, and the reader is referred to 
it (and references therein) for further details. In short, 
Hipparcos/Tycho photometry is adopted as a fair rep- 
resentation of typical Gaia photometric data, and de- 
voted ground-based observations obtained with the Asiago 
1.82m telescope are used to simulate expected Gaia spec- 
troscopic data. The latter are obtained in the same 
Gaia wavelength region (AA8480— 8740 A), at a higher 
resolution (R=20000) than is planned for the satel- 



lite (R=11500) to compensate for the lower number of 
recorded radial velocity points (~30 vs. the ~100 expected 



Table 2. Summary of Hipparcos {Hp) and Tycho (Bt, 
Vt) photometric data and ground based radial velocity 
observations: number of data, mean S/N and mean stan- 
dard error (magnitudes for photometry, km sec~^ for ra- 
dial velocities) for each of the three program stars. 





Hip 




Tycho 




RV 






N a{Hp) 


N 


(j{Bt) o{Vt) 


N 


S/N cr 


(RV) 


SV Cam 


113 0.021 


167 


0.34 0.31 


35 


40 


10 


BS Dra 


88 0.016 


143 


0.20 0.19 


27 


36 


6 


HP Dra 


81 0.012 


122 


0.11 0.09 


29 


62 
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from Gaia). Gaia data will obviously largely surpass the 
Hipparcos/Tycho photometry adopted here in terms of 
number of epochs, number and diagnostic capabilities of 
the photometric bands, and therefore even better perfor- 
mance than we have found here can be expected from the 
satellite. 

2. The targets 

As in the previous papers in the series (Munari et al. 2001, 
Paper I, Zwitter et al. 2003, Paper II and Marrese et al. 
2004, Paper 111), we have selected both newly discovered 
eclipsing binaries (i.e., those lacking spectroscopic and 
photometric orbit solutions in the literature), as well as bi- 
naries with already published orbital solutions (although 
not with data in the Gaia spectral range) to permit a com- 
parison. In our analysis, to properly simulate Gaia, we 
have ignored all pre-existing ground-based data, our solu- 
tions resting exclusively on Hipparcos/Tycho photometry 
and Gaia spectroscopy as collected in Asiago. Similarly, we 
have intentionally avoided selecting initial parameters for 
our modeling trials from published solutions. The three 
targets for this fourth paper in the series are SV Cam, 
BS Dra and HP Dra. The first one is still a controversial 
object in spite of hundreds of papers devoted to it, while 
some orbital modeling has been already attempted for the 
second. The third one is a relatively unstudied eclipsing 
system. Some basic quantities are given in Table 1. 

SV Cam. This is reported in literature as a mid- 
F to mid-G system with 0.6 day period and a light- 
curve of the (3 Lyrae type. There is a long history of 
observations and analyses of this system, with extensive 
photometry (e.g., Kjurkchicva, Marchcv & Ogloza 2000; 
Patkos 1982; Zeilik ct al. 1988), and spectroscopy (e.g., 
Hempelmann et al. 1997; Kjurkchieva, Marchev & Zola 
2002; Lehmann, Hempelmann & Wolter 2002; Ozeren et 
al. 2001; Pojmaiiski 1998; Popper 1996; Rainger, Hilditch 
& Edwin 1991 and Rucinski ct al. 2002). Most authors 
have interpreted asymmetries in the light curves and the 
variation in them as due to spots. With the inclusion 
in models of high latitude, long-lasting spots, it is not 
surprising that there is substantial disagreement about 
the spectral type of the hotter component, with esti- 
mates ranging from F5 to G8 (Popper 1996). The dif- 
ference in cK'lipse depths shows that the secondary star, 
as we refer to the component eclipsed at secondary min- 
imum, is much fainter than the primary (cf. Figure 2) 
and therefore it is not a surprise that it is not easily de- 
tectable and measurable in our spectra. A similar conclu- 
sion was reached by Popper (1996), based on the photome- 
try of Zeilik et al. (1988) and by Patkos and Hempelmann 
(1994). Consequently we include in the analysis only the 
radial velocity curve of the hotter, more luminous com- 
ponent. This means that, from our data, the mass ratio 
is not independently determined by spectroscopic means, 
and its determination must depend on the curvature of 
the maxima in the light curves, which are relatively noisy. 
This system thus becomes an interesting test case for Gaia 
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analyses, both because of activity and because of the low 
luminosity of the secondary. 

BS Dra. This system is a partially eclipsing Algol sys- 
tem with similar components and a period of ~3.4 days. 
The Gaia-like spectra reveal two stars of nearly identical 
spectral types and luminosities. A comparison with the 
spectral atlases by Munari & Tomasella (1999), Marrese 
et al. (2003) and Zwitter et al. (2004) obtained in the 
same wavelength region at the same (R=20 000) resolu- 
tion, indicates that the the spectral type is ~ F3V and 
a metallicity lower than solar ([Fe/H]~ —0.4). Previous 
studies were carried out by Popper (1971), Giidiir et al. 
(1979) and by Russo et al. (1981). 

HP Dra. This system was discovered as an eclipsing 
binary by Hipparcos, but the reported period (6.67 days) 
is incorrect. We find from our spectroscopic data that this 
Algolid system has a period of ~ 10.761 days. A similar 
value has been found also by Kurpinska-Winiarska et al. 
(2000). A comparison with the above cited spectral atlases 
suggests a spectral type of ~ F9V and a solar metallicity. 
The lack of data at minimum light in the Hipparcos and 
Tycho data sets proved to be a severe problem in modeling 
this system, as discussed below. 

3. The data 

3.1. Spectroscopic observation and radial velocities 

The Gaia simulating spectroscopic observations were ob- 
tained with the Echelle and CCD spectrograph on the 

I. 82 m telescope operated by Osservatorio Astronomico 
di Padova atop Mt. Ekar (Asiago), the same facility and 
instrument as used for the targets of the previous papers 
in this series. The same wavelength range (8480-8740 A), 
reciprocal linear dispersion (0.25 A/pix), and resolution 
(R = A/AA = 20 000) as for the data reported in Papers I, 

II, and III similarly apply. F\irther details of the spectro- 
graph and the spectroscopy can be found in these earlier 
papers. The diary of observations of the stars discussed in 
this paper appear in Table 3. Sample spectra in the Ca II 
triplet region for each of the three systems are shown in 
Fig. 1. The number of spectroscopic observations and their 
accuracy are summarized in Table 2. 

3.2. Hipparcos/Tycho photometry 

Hipparcos and Tycho epoch photometry was retrieved 
from CDS. The conversion of Hipparcos and Tycho photo- 
metric quantities to the Johnson-Cousins system, follows, 
as in previous papers, the transformation equations pro- 
vided in the Hipparcos Catalogue (ESA 1997): 

Vj ^ Vt - 0.090 X {B-V)t, (1) 

{B-V)j = 0.850 X (B-V)t (2) 

Details on the number of observations and their accuracy 
are given in Table 2. The relationships among {B — V), 
spectral type, temperature, and bolometric correction 
were taken from Popper (1980). 
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Table 3. Journal of radial velocity data. The columns give the Heliocentric Julian Day and the radial velocities in 
km sec~^ for each of the three systems discussed in this paper. 



4. The modeling software 

We began by using the Wilson-Devinney base program 
(Wilson 1998), with modifications noted by Kallrath et 
al. (1998) and Kalhath & Milone (1999). Since Paper I, 
the University of Calgary version has further evolved to 
wd98k93h, a Unix version which not only makes use of 
Kurucz stellar atmosphere models applied to each grid el- 
ement, but also has sufficiently numerous grid elements 
that transit eclipses by objects as small as a Saturn-sized 
planet can be modeled. This version of wd98k93 has been 
incorporated into the WD2002 package that runs under 
Linux or Microsoft Windows. The WD2002 package per- 



mits self-iterated, damped least squares differential cor- 
rections, as did its predecessors, and now permits looping 
over a range of a single parameter, such as q or i, to im- 
prove optimization in highly correlated cases. 

The program requires tables of the ratio of Kurucz at- 
mospheres to black-body fluxes in each of the Hp, Bt 
and Vt passbands as well as in a square passband cen- 
tered on the Ca II triplet region, used in connection with 
the radial velocity curves. In addition the limb-darkening 
coefficients for these passbands are needed. The former 
were provided by C.R. Stagg and MDW, and the latter 
by W. Van Hamme together with the correct limb dark- 
ening coefficients for appropriate effective temperatures. 
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surface gravities and metallicities. Additionally a two-way 
coefficient interpolator was provided by D. Terrell. 

The usual operation begins with reasonable guesses for 
the initial parameters based primarily on the spectrosc- 
copic data and secondarily on the colour indices of the 
Tycho observations. With WD2002, we try to ascertain 
the correct global minimum region of parameter space, 
by means of simplex iterations which explore all the pa- 
rameter space. Then final parameters are found by iter- 
ated damped least-squares DC runs. When convergence is 
achieved, the process is concluded. The uncertainties in 
the orbital parameters quoted in Table 4 come from final 
modeling check runs carried out with wd98k93h. 

5. Analyses 

5.1. Procedure 

All modeling runs were carried out with time-ordered 
data, so that period and epochs were adjustable param- 
eters. Initial periods and epochs were determined from 
the spectroscopy, as were the spectral types of the com- 
ponents (except for the undetectable secondary star of 
SV Cam, as we note in the previous section). Typically, 
the following parameters were adjusted: semi-major axis 
(a), systemic radial velocity {V^), inclination (i), temper- 
ature of the secondary star, the star at inferior conjunc- 
tion at the designated primary minimum, (T2), the mod- 
ified Kopal potentials of each star {^1^2), the mass ratio 
{q — M2/M1), the passband luminosity in units of An (Li) 
and, as already mentioned, the epoch (To) and period (P). 
Additional suites of parameter adjustments were carried 
out for BS Dra and HP Dra. Convective atmosphere val- 
ues for gravity brightening and albedo coefficients g = 0.32 
and A = 0.500, respectively, were assumed for stars sim- 
ilar to or cooler than the Sun, and radiative atmosphere 
values of 1 for both g and A were assumed for higher tem- 
perature stars. A full-precision grid number of 30 was used 
for all stars. 

5.2. Solutions 

5.2.1. SV Camelopardalis 

The temperature adopted for the hotter star, Ti = 5848 K, 
was based on a compromise. The spectra in the Ca II 
triplet region indicates a temperature ^ 5800 K or cooler. 
The Tycho colors, however, support an slightly earlier 
spectral type. The star is rapidly rotating, evolved and 
sports numerous surface spots (cf Kjurkchieva et al. 2002 
and references therein) that are responsible for a signifi- 
cant uncertainty when comparing with template field stars 
that usually do not have spotty surfaces and rotate more 
slowly. The adopted 5848 K corresponds to a spectral type 
of G2-G3, as suggested also by Rainger et al. (1991). The 
temperature of the secondary component is determined 
solely from the light curve analysis. 

The single radial velocity curve and the Hipparcos light 
curve provide the bulk of the modeling input informa- 
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0.0 0.5 1,0 phase 



Fig. 2. Hipparcos Hp and Tycho Vr, Bt, {B - V)t 
lightcurves and the radial velocity curve of SV Cam folded 
onto a period P ^ 0.593 days. The lines represent the so- 
lution given in Table 4. 



tion because of large scatter in the Tycho light curves. 
The analysis was performed fully on two different mod- 
els. The first model relies on only the radial velocity and 
-ffp measurements, while the second includes in addition 
Tycho Bt and Vt data. Although the formal fitting er- 
ror is smaller for the first model, the determination of the 
mass ratio in this case depends strongly on the curvature 
of the light curve outside of eclipse; therefore, notwith- 
standing the greater scatter in the Tycho data, three light 
curves may provide a better determination of this critical 
quantity than one, with appropriate curve weights relating 
to the intrinsic scatter in the light curves. 
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Table 4. Modeling solutions. The uncertainties are formal mean standard errors of the solution. The last four rows give 
the r.m.s of the observed points from the derived orbital solution. The estimation of error in the adopted temperature 
for the primaries is ~150 K. 



parameter (units) 


bV Cam 


rSb Ura 


rlF Ura 


Period (days) 


U.oyaU7oa7±U.UUUUUUoD 


im c 1 J n An AAA CI 
o.oD4(Ji54±U. UUUUUoi 


1A Ti^lC/l-LA AAAAA 


iijpocn [tiJU) 


O/i/iTinT Qno/i_Ln aai o 


O/I/IOKAO QAQQC-LA AAAQA 


0/1/1TA/10 OC/1_LA AAQ 

z44 / y4z.ZD4±U.UUo 




ci.41±0.1ci 


12.94±0.25 


2d.4o±0.43 


(km sec J 


— 9.2±1.4 


— 0.5±1.5 


1 A"? 1 A 1^ 

— 1d.07±0.d3 




0.750±0.053 


0.986±0.038 


0.988±0.028 


i (deg) 


73.8±1.4 


89.54±0.14 


87.87±0.14 


e 


0.0 


0.0 


0.057±0.010 


(radians) 






1.14±0.15 


du)/dt (radians day~^) 






-0.000156±0.000080 


Ti (K) 


5848 


6618 


6000 


T2 (K) 


4061±176 


6626±153 


6386±207 


A(Ti-T2) (K) 


1787±93 


-8±25 


-386±143 


Hi 


4.29±0.33 


10.14±0.14 


23.81±0.65 


^2 


3.39±0.10 


10.15±0.46 


31.1±2.6 


»1 (R0) 


0.98±0.10 


1.415±0.022 


1.172±0.034 


»2 (R©) 


1.18±0.12 


1.396±0.071 


0.887±0.086 


Ml (Mq) 


0.862±0.099 


1.294±0.079 


1.099±0.023 


il J /TV /f N 

M2 (M0) 


0.d4o±0.087 


1.27D±0.0y2 


1.102±0.02d 




4.69±0.25 


3.34±0.10 


4.18±0.13 




5.86±0.29 


3.37±0.15 


4.51±0.24 


log g\ (cgs) 


4.40±0.14 


4.248±0.040 


4.341±0.034 


log 91 (cgs) 


4.11±0.15 


4.254±0.075 


4.584±0.086 


cTi (avg.wt.) 


±0.0545 


±0.0109 


±0.0130 


o'w,i,2 (km sec~^) 


±0.66 


±0.19, 0.15 


±0.17, 0.21 


(mag) 


±1.25 






(TVt (mag) 


±0.79 






(JHp (mag) 


±0.26 


±0.17 


±0.09 



The absolute parameters and other derived quantities 
for the adopted second model are given in Table 4 and the 
fittings are shown in Figure 2. SV Cam is indeed a chal- 
lenging single-lined spectroscopic binary, nevertheless we 
obtained a reasonable solution with accuracies of ~13% 
in the masses and ~10% in the radii. The cooler star con- 
tributes little to B and V light curves and even in the 
Ca II triplet region contributes less than 15%, according 
to our model, which indicates that it is very close to its 
Roche lobe. 

The photometric solution of Albayrak ct al. (2001) 
rests on a temperature of 6440 K for the hotter com- 
ponent and sets of 242 observations each for BVR pass- 
bands. While their period and orbital separation are quite 
close to oins in Table 4 (they quote _P=0. 593071 days and 
a=3.87±0.07 i?©), their mass ratio ((7=0.56) and inclina- 
tion (i=89.6±0.8) are quite inconsistent with ours. There 
are many possible reasons for the differences: inclusion or 
exclusion of RV data, different light curve data, different 
limb-darkening coefficients, and the need for treatment of 
asymmetries. The Hipparcos/Tycho light curves present 
no explicit evidence of asymmetries that would justify spot 
modeling, which was therefore not included in our analy- 



ses, highlighting a possible shortcoming in Gaia studies of 
active surface, single lined eclipsing binaries. 

Our assessments of the SV Cam system are based ex- 
clusively on the treatment of our adopted input data. As 
our anonymous referee has pointed out, these assessments 
arc not necessarily correct. Indeed, higher quality ground- 
based data on SV Cam and higher resolution spectroscopy 
has already indicated that the secondary component has 
been observed and that the mass ratio is very likely 0.6 
(somewhat different from our determination, described be- 
low). Although neither the precision in the Hipparcos- 
Tycho photometry nor the spectral resolution in our RV 
data match those of other studies, our modeling tools per- 
mit distortion in both curves to be treated adequately, 
given any discernible spot regions. Therefore our failure 
to find spot regions may be an indication of a relatively 
"quiet" interval in the system's behaviour, when the ac- 
tivity was below the detection level of the photometry. In 
such mismatch may arise because of differences 

in activity level at the separate epochs at which the data 
were acquired. The variability impacts our conclusions, 
and divergcncx; from previous ground-based investigations 
indicates that the results for SV Cam are very important 
for GAIA mission planning. 
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Fig. 3. Hipparcos Hp and Tycho Vr, Bt, {B - V)t: 
lightcurves of BS Dra folded onto a period P = 3.364 
days. The lines represent the solution given in Table 4. 



ponents. Moreover, the additional information about the 
curvature at the light curve maxima was not important 
to the determination of the mass ratio, imlikc the case for 
SV Cam. It was found that the logarithmic form of the 
limb- darkening appeared to improve the fits slightly, and 
these were retained for the analyses. 

There are very few Hipparcos/Tycho photometric 
points covering the eclipse phases. One of the eclipses 
in the H-p light curve has no data points within 10% of 
minimum value and the other has only three points alto- 
gether, all of which are on the descending branch. Yet, 
if we assume symmetry, the shape and the depths are 
defined, even if not with great precision. The accuracy 
of the Hipparcos data for this star (uffp =0.016) is high 
enough to permit the few points in the light curve minima 
to be highly weighted to insure that the computed light 
curve fits them closely. All other points were individually 
weighted inversely to their standard errors, and scaled. 
The determination of the modified Kopal potentials and 
thus the radii depend critically on this close-fitting pro- 
cess. 

Both the mass ratio and temperature ratios are near 
1. During modeling, both the mass ratio and the tem- 
perature ratio switched repeatedly across this boundary. 
Solutions with either star as the primary component pro- 
duce nearly the same fit error, but the adopted converged 
solution has the epoch given in the Hipparcos catalogue. 
This epoch is, however, half a cycle different from that 
of Popper (1971). The alternate model, where the compo- 
nents are interchanged, converged only with the assump- 
tion of a small eccentricity, but under this circumstance, 
the argument of periastron, oj, could not be determined. 
The stellar and system parameters arc given in Table 4 
for the adopted model with [Fe/H] = —0.5. Our results 
can be compared with those of Russo et al. (1981), who 
used an earlier version of the Wilson-Devinney program 
to analyze previously published BV observations, proba- 
bly those of Giidiir et al. (1979). Their values for period, 
orbital separation, inclination and modified Kopal poten- 
tials are close (even if not always within Icr level) to our 
results. 



5.2.2. BS Draconis 

The inspection of the spectra of BS Dra suggests an F3 
spectral type and a metallicity lower than solar. The tem- 
perature of the primary was fixed to 6619 K (Popper 1980) 
and the modeling was carried out with [Fe/H]=0.0, —0.5, 
and —1.0 metallicities. Although scarcely significant, the 
fitting was seen to be slightly better for the —0.5 metal- 
licity, as expected for its proximity with the [Fe/H]=— 0.4 
estimated from our spectra. 

The Tycho and Vt hght curves are so noisy that 
they essentially add nothing to the modeling, and were 
excluded from the analysis. Dropping the information on 
color variation does not affect the accuracy of the solution 
given the almost identical temperature for the two com- 



5.2.3. HP Draconis 

HP Dra revealed itself as the most challenging of the 12 ob- 
jects in this series of papers. The HP Dra components ap- 
pear to be similar in temperature and luminosity. Because 
of their great scatter and sparseness at minimum light, 
we rcinioved the Tycho data from the modeling sets, even 
though the minima are weakly discernible, and modeled 
the Hipparcos and radial velocity data alone. There are 
even fewer data points in the minima for HP Dra than 
there are for BS Dra (9 vs. 11). The placements of the 
data points at minimum light are critical to the determi- 
nation of the radii and potentials of the stars. Although 
the precision of the Hipparcos data is good (cthp =0.012), 
it still proved very difficult to achieve a final, fully conver- 
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Fig. 4. Hipparcos Hp and Tycho Vt, Bt, (B — V)t 
lightcurves of HP Dra folded onto a period P = 10.761 
days. The lines represent the solution given in Table 4. 



gent solution. The major trouble is that the non-coeval 
light and radial velocity curves (about 10 years elapsed 
between Hipparcos and Asiago observations) appear to be 
not fully compatible. If this problem is indeed connected to 
the non-coeval data sets, this would not be a problem for 
Gaia where all astrometric, photometric and spectroscopic 
observations will strictly share the same epochs. Several 
strategies were attempted to deal with this circumstance. 

1. Initially, the radial velocities were modeled indepen- 
dently of the photometry, with assumed temperatures 
of ~ 6000 K for each component, appropriate for the 
spectral type (~ F9) and with equal luminosities. This 
yielded elements that produced an optimal fit to the 
radial velocity data, although the assumptions about 



some of the adopted parameters (such as ^^1,2) could 
not be rigorously checked. The RV modeling indicated 
the need to adjust the eccentricity and argument of 
periastron. A slight difference in mass between the 
components was also found. The parameters which 
were adjusted in these trials were: e, lu, V-y, q, Tq, 
and P, which are, respectively, the eccentricity, argu- 
ment of periastron, system radial velocity, mass ratio 
(M2/M1), epoch, and period. 

2. Further attempts at modeling the combined light and 
RV curves with these parameters fixed were also at- 
tempted. The results gave minima which were system- 
atically too shallow. Attempts to weight points in the 
minima more heavily also failed to produce conver- 
gence. 

3. Adjustment of the third light, which failed to improve 
the fitting errors. 

4. The phased Hp light curve data appear to be mod- 
ulated by a rough sinusoid at maximum light. This 
effect is also seen in the 1989-1990 V light curve of 
Kurpinska-Winiarska et al. (2000). Their light curve 
shows relatively little scatter when phased with a sim- 
ilar but independently determined period. One possi- 
bility is a reflection effect of some sort. This is not 
likely, given the scale of the orbit and similarity of 
the temperatures of the component stars, but it was 
investigated anyway. Introduction of a more compli- 
cated reflection effect option in the Wilson-Devinney 
program also failed to provide adequate fittings to the 
Hp light curve, even when the size of the primary 
star was initially increased (by decreasing its poten- 
tial). The computing time increases greatly when this 
option is adopted in treating an eccentric orbit case, 
such as HP Dra, so we performed only one trial with 
the enhanced (2-pass) reflection effect. The fitting was 
not quite as good as for the simple reflection case, in 
which the illuminating star is assumed to be a point 
source and only one reflection is treated. 

5. Because intrinsic variation may also be affecting the 
light curve, we attempted to introduce a large starspot 
onto the visible hemisphere of the component at in- 
ferior conjunction during eclipses: either a hot spot 
on the 'secondary' star visible at primary minimum 
(phase 0.0) or a cool, spot on the 'primary' star vis- 
ible at secondary minimum (phase 0.5). With spots 
centered on a longitude of 180° and with a radius of 
90°, the modulation was adequately fit but a solution 
with this artifice still failed to converge, in either spot 
placement case. The existence of a spot on stars just 
slightly hotter than the Sun would not be surprising, 
but we have not found lower fit errors by assuming 
convective envelopes in these components (although 
not every model was tested), thus star spots appear 
to be both unnecessary and insufficient to model this 
system, given only the current suite of data that have 
been modeled. 

6. An investigation of variability in the residuals failed to 
show periodicities differing slightly from the derived 
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eclipsing period, viz., over the range 3.364000000 to 
3.364033141 days. 

7. Adjustment of dP/dt equally failed to produce an im- 
provement. 

8. On the contrary, adjustment of the dw/dt term did 
produce an improvement in fitting of both light curve 
and radial velocity curves (when modeled simultane- 
ously). This finally yielded a converged solution that 
allowed us to combine into a coherent picture both 
photometric data and radial velocities in spite of the 
10 year elapsed between them. This is the adopted so- 
lution presented in Table 4. 

Even though the system seems to be at or just be- 
yond the data-paucity limit for analysis of partially eclips- 
ing systems, we were able to achieve a fully convergent 
and self-consistent solution. It is worth repeating here 
that the two major problems encountered in the analy- 
sis of HP Dra (i.e., non-coevality of photometric and ra- 
dial velocity curves and paucity of the photometric points 
covering the eclipse phases) will not pertain to Gaia (cf. 
Katz et al. 2004). The satellite will obtain simultaneous 
spectroscopic and photometric data during the planned 
5 year mission, the photometry at about ~355 indepen- 
dent epochs (cf. Munari et al. 2004) compared to the ~100 
Hipparcos epochs. 

Nevertheless, the resulting model cannot be said to be 
a truly accurate one. With only three points in one of 
the minima (and only one point near minimum value), it 
is clear that the relative surface brightnesses of the two 
stars cannot be obtained from the Hipparchos photome- 
try alone, the depth providing less information than the 
branches of the minima, which the model must fit. In fact 
the light curve produced by Kurpinska-Winiarska et al. 
(2000) shows the secondary minimum to be the shallower. 
This suggests that there is indeed a temperature differ- 
ence between the components, but in the opposite sense 
from that found here. The large uncertainty in the tem- 
perature difference is therefore reasonable. The difference 
in the potentials is determined largely by the phases of 
contact and is therefore also uncertain, but the outer con- 
tacts are not badly defined by the data. Our model indi- 

Table 5. Comparison between the Hipparcos distances 
and those derived from the parameters of the modeling 

solution in Table 4 assuming an uncertainty of ±150 K in 
the temperature adopted for the primary in each system. 

Hipparcos this paper 





(pc) 


(pc) 


SV Cam 


85?i 


87 ±8 


BS Dra 


208?|f 


172 ±g 


HP Dra 


80fi 


73 ±4 



cates that there is a small orbital eccentricity, in rough 
agreement with a preliminary value (e=0.043) found by 
Kurpinska-Winiarska et al. (2000). As far as we know, 
there has been no evidence presented for a third body 
in the system, but we could not achieve convergence in 
the modeled data set without adjustments of the apsidal 
motion. For HP Dra, complete light- and radial velocity 
curves exist and the full analysis promised by Kurpinska- 
Winiarska et al. (2000) is awaited with great interest. If 
the apsidal motion mentioned here can be confirmed, ad- 
ditional photometry across a range of bandpasses, times 
of minima, and spectroscopy, should be pursued. 

It is possible that further spot modeling on a fuller data 
set will be able to resolve the question of the true nature 
of HP Dra. For now, we indicate only the possibility of 
apsidal motion, and the limitations for GAIA models given 
our current suite of tools to test that mission's capability. 

6. Conclusions 

The distances derived from the orbital solutions are com- 
pared to those from the Hipparcos Catalogue in Table 5. 
Despite the fairly large uncertainties in the fundamen- 
tal parameters, the derived distances to all three sys- 
tems agree with the astromctric distances from Hipparcos- 
Tycho within or very close to their one-sigma errors. 

The overall agreement confirms similar findings in pre- 
vious papers in this series and together they reinforce con- 
fidence in the overall quality and astrophysical potential 
of orbital modeling of eclipsing binary data to be obtained 
by Gaia. 

As in other papers in this series, we have been discrim- 
inate in the use of the Tycho mission data, because of its 
limited precision, and effectively used them in the present 
paper only for the SV Cam system for reasons mentioned 
earlier. For fainter systems, the Tycho data are useful only 
for the limited mean color information that they provide. 
For the BS Dra and HP Dra systems, they proved unsuit- 
able for modeling purposes. 

Although we conclude that the results for the three 
systems presented here are plausible, in context with ex- 
isting, higher quality data other than our adopted data 
suite, we note the possibility of systematic error in the 
analyses based on photometry only from Hipparcos and 
Tycho missions, even when these are coupled with spec- 
troscopy that closely matches GAIA's. In particular, the 
results for SV Cam, and possibly also for HP Dra, point 
to the need to diagnose spot regions in systems detected 
by GAIA. Of course, it is possible that a system will ap- 
pear " quiet" during times of data acquisition, but GAIA 
will acquire all its types of data at the same instant, so 
they will be self-consistent, at least. Moreover, the suite of 
photometric passbands recently selected for GAIA permits 
the discrimination of system properties from photometric 
indices of passbands across a broad region of the spec- 
trum, making the color effects of spotted regions easier to 
detect. 
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